Interindividual behavioural variation in response to elevated CO2 predicts mRNA transcript abundance of genes related to acid-base regulation in medaka (Oryzias latipes).

Rising carbon dioxide (CO2) in aquatic ecosystems due to climate change is a challenge for aquatic ectotherms. We examined whether interindividual variation in behavioural responses to CO2 could predict how a teleost fish would respond to elevated CO2 for multiple phenotypic and molecular traits. To this end, we first quantified behavioural responses of individuals exposed to acute elevated CO2, and used these to assign individuals as either high or low responders relative to the population mean. Subsequently, we exposed both high and low responders to elevated CO2 for 6 weeks and quantified the effect on body condition, behaviour, and mRNA transcript responses of gill and liver genes associated with relevant physiological processes. Generally, we found few relationships between the phenotypic groups and body condition and behaviour following the CO2 exposure period; however, stark differences between the phenotypic groups with respect to gene transcripts from each tissue related to various processes were found, mostly independent of CO2 exposure. The most pronounced changes were in the gill transcripts related to acid-base regulation, suggesting that the observed behavioural variation used to assign fish to phenotypic groups may have an underlying molecular origin. Should the link between behaviour and gene transcripts be shown to have a fitness advantage and be maintained across generations, interindividual variation in behavioural responses to acute CO2 exposure may be a viable and non-invasive tool to predict future population responses to elevated aquatic CO2.

A high-quality chromosome-level genome assembly of the Chinese medaka Oryzias sinensis.

Oryzias sinensis, also known as Chinese medaka or Chinese ricefish, is a commonly used animal model for aquatic environmental assessment in the wild as well as gene function validation or toxicology research in the lab. Here, a high-quality chromosome-level genome assembly of O. sinensis was generated using single-tube long fragment read (stLFR) reads, Nanopore long-reads, and Hi-C sequencing data. The genome is 796.58 Mb, and a total of 712.17 Mb of the assembled sequences were anchored to 23 pseudo-chromosomes. A final set of 22,461 genes were annotated, with 98.67% being functionally annotated. The Benchmarking Universal Single-Copy Orthologs (BUSCO) benchmark of genome assembly and gene annotation reached 95.1% (93.3% single-copy) and 94.6% (91.7% single-copy), respectively. Furthermore, we also use ATAC-seq to uncover chromosome transposase-accessibility as well as related genome area function enrichment for Oryzias sinensis. This study offers a new improved foundation for future genomics research in Chinese medaka.

Sexually dimorphic dynamics of the microtubule network in medaka (Oryzias latipes) germ cells.

Gametogenesis is the process through which germ cells differentiate into sexually dimorphic gametes, eggs and sperm. In the teleost fish medaka (Oryzias latipes), a germ cell-intrinsic sex determinant, foxl3, triggers germline feminization by activating two genetic pathways that regulate folliculogenesis and meiosis. Here, we identified a pathway involving a dome-shaped microtubule structure that may be the basis of oocyte polarity. This structure was first established in primordial germ cells in both sexes, but was maintained only during oogenesis and was destabilized in differentiating spermatogonia under the influence of Sertoli cells expressing dmrt1. Although foxl3 was dispensable for this pathway, dazl was involved in the persistence of the microtubule dome at the time of gonocyte development. In addition, disruption of the microtubule dome caused dispersal of bucky ball RNA, suggesting the structure may be prerequisite for the Balbiani body. Collectively, the present findings provide mechanistic insight into the establishment of sex-specific polarity through the formation of a microtubule structure in germ cells, as well as clarifying the genetic pathways implementing oocyte-specific characteristics.

Evidence for hybridization-driven heteroplasmy maintained across generations in a ricefish endemic to a Wallacean ancient lake.

Heteroplasmy, the presence of multiple mitochondrial DNA (mtDNA) haplotypes within cells of an individual, is caused by mutation or paternal leakage. However, heteroplasmy is usually resolved to homoplasmy within a few generations because of germ-line bottlenecks; therefore, instances of heteroplasmy are limited in nature. Here, we report heteroplasmy in the ricefish species Oryzias matanensis, endemic to Lake Matano, an ancient lake in Sulawesi Island, in which one individual was known to have many heterozygous sites in the mitochondrial NADH dehydrogenase subunit 2 (ND2) gene. In this study, we cloned the ND2 gene for some additional individuals with heterozygous sites and demonstrated that they are truly heteroplasmic. Phylogenetic analysis revealed that the extra haplotype within the heteroplasmic O. matanensis individuals clustered with haplotypes of O. marmoratus, a congeneric species inhabiting adjacent lakes. This indicated that the heteroplasmy originated from paternal leakage due to interspecific hybridization. The extra haplotype was unique and contained two non-synonymous substitutions. These findings demonstrate that this hybridization-driven heteroplasmy was maintained across generations for a long time to the extent that the extra mitochondria evolved within the new host.

Deciphering the origin of developmental stability: The role of intracellular expression variability in evolutionary conservation.

Progress in evolutionary developmental biology (evo-devo) has deepened our understanding of how intrinsic properties of embryogenesis, along with natural selection and population genetics, shape phenotypic diversity. A focal point of recent empirical and theoretical research is the idea that highly developmentally stable phenotypes are more conserved in evolution. Previously, we demonstrated that in Japanese medaka (Oryzias latipes), embryonic stages and genes with high stability, estimated through whole-embryo RNA-seq, are highly conserved in subsequent generations. However, the precise origin of the stability of gene expression levels evaluated at the whole-embryo level remained unclear. Such stability could be attributed to two distinct sources: stable intracellular expression levels or spatially stable expression patterns. Here we demonstrate that stability observed in whole-embryo RNA-seq can be attributed to stability at the cellular level (low variability in gene expression at the cellular levels). We quantified the intercellular variations in expression levels and spatial gene expression patterns for seven key genes involved in patterning dorsoventral and rostrocaudal regions during early development in medaka. We evaluated intracellular variability by counting transcripts and found its significant correlation with variation observed in whole-embryo RNA-seq data. Conversely, variation in spatial gene expression patterns, assessed through intraindividual left-right asymmetry, showed no correlation. Given the previously reported correlation between stability and conservation of expression levels throughout embryogenesis, our findings suggest a potential general trend: the stability or instability of developmental systems-and the consequent evolutionary diversity-may be primarily anchored in intrinsic fundamental elements such as the variability of intracellular states.

Ran-GTP assembles a specialized spindle structure for accurate chromosome segregation in medaka early embryos.

Despite drastic cellular changes during cleavage, a mitotic spindle assembles in each blastomere to accurately segregate duplicated chromosomes. Mechanisms of mitotic spindle assembly have been extensively studied using small somatic cells. However, mechanisms of spindle assembly in large vertebrate embryos remain little understood. Here, we establish functional assay systems in medaka (Oryzias latipes) embryos by combining CRISPR knock-in with auxin-inducible degron technology. Live imaging reveals several unexpected features of microtubule organization and centrosome positioning that achieve rapid, accurate cleavage. Importantly, Ran-GTP assembles a dense microtubule network at the metaphase spindle center that is essential for chromosome segregation in early embryos. This unique spindle structure is remodeled into a typical short, somatic-like spindle after blastula stages, when Ran-GTP becomes dispensable for chromosome segregation. We propose that despite the presence of centrosomes, the chromosome-derived Ran-GTP pathway has essential roles in functional spindle assembly in large, rapidly dividing vertebrate early embryos, similar to acentrosomal spindle assembly in oocytes.

Loss of function in somatostatin receptor 5 has no impact on the growth of medaka fish due to compensation by the other paralogs.

Somatic growth in vertebrates is regulated endocrinologically by the somatotropic axis, headed by the growth hormone (GH) and the insulin growth factor-I (IGF-I). Somatostatin (Sst), a peptide hormone synthesized in the hypothalamus, modulates GH actions through its receptors (Sstr). Four Sstr subtypes (Sstr 1-3 and 5) have been identified in teleosts. However, little is known about whether they have a specific function or tissue expression. The aim of this study was to determine the role of sstr2 and sstr5 in the growth of the medaka (Oryzias latipes). The assessed expression pattern across diverse tissues highlighted greater prevalence of sstr1 and sstr3 in brain, intestine and muscle than in pituitary or liver. The expression of sstr2 was high in all the tissues tested, while sstr5 was predominantly expressed in the pituitary gland. A CRISPR/Cas9 sstr5 mutant with loss of function (sstr5-/-) was produced. Assessment of sstr5-/- indicated no significant difference with the wild type regarding growth parameters such as standard length, body depth, or peduncle depth. Furthermore, the functional loss of sstr5 had no impact on the response to a nutritional challenge. The fact that several sstr subtypes were upregulated in different tissues in sstr5-/- medaka suggests that in the mutant fish, there may be a compensatory effect on the different tissues, predominantly by sstr1 in the liver, brain and pituitary, with sstr2 being upregulated in pituitary and liver, and sstr3 only presenting differential expression in the brain. Analysis of the sstr subtype and the sstr5-/- fish showed that sstr5 was not the only somatostatin receptor responsible for Sst-mediated Gh regulation.

Comprehensive analysis of diel rhythmic expression of the medaka toll-like receptor gene family.

Several immune-related genes, including Toll-like receptors (TLR), are associated with circadian rhythms in mammals. However, information on the circadian rhythmic expression of TLRs in fish is limited. In this study, we aimed to analyze the regulation of diel oscillations in the expression of TLR genes in Japanese medaka (Oryzias latipes). The expression analysis revealed diel expression patterns of tlr1, tlr5m, tlr21, and clock genes (bmal1 and clock1) under a 12 h light:12 h dark cycle. The clock gene response element (E-box) was identified in the transcriptional regulatory regions of tlr1, tlr5m, and tlr21. Moreover, overexpressed bmal1 and clock1 enhanced expression levels of tlr1, tlr5m, and tlr21 in medaka embryo (OLHdrR-e3) cells. The expression of tlr1, tlr5m, and tlr21 was significantly decreased in OLHdrR-e3 after generating a bmal1 knockdown using a morpholino oligo. These results indicate the regulation of the diel rhythmic expression of several fish TLRs by clock genes.

Effect of diisobutyl adipate on the expression of biomarker genes that respond to endocrine disruption and on gonadal sexual differentiation in Japanese medaka (Oryzias latipes).

Phthalate and non-phthalate plasticizers are used in polymer materials, such as plastic and rubber. It has recently been found that diisobutyl adipate (DIBA), which is considered an environmentally safe non-phthalate plasticizer, potentially acts as a thyroid disruptor in fish. Here, we investigated the sexual hormone effects of DIBA based on the expression levels of genes that respond to endocrine disruption and sexual hormone activity in the livers and gonads, and on gonadal sexual differentiation in Japanese medaka. Compared with the control group, the mRNA expression of chgH, vtg1, vtg2, and esr1 was significantly suppressed in the livers of DIBA exposed XX individuals. Furthermore, the mRNA expression of gsdf was significantly upregulated and downregulated in the gonads of XX and XY individuals, respectively. The mRNA expressions of esr1 and esr2b were significantly suppressed by DIBA exposure in the gonads of both XX and XY individuals. These observations suggest that DIBA has potential androgenic activity in Japanese medaka. However, normal testes and ovaries were observed in respective XY and XX medaka after DIBA exposure; therefore, these results suggest that DIBA may have weak androgenic activity.

A transcriptional program underlying the circannual rhythms of gonadal development in medaka.

To cope with seasonal environmental changes, organisms have evolved approximately 1-y endogenous circannual clocks. These circannual clocks regulate various physiological properties and behaviors such as reproduction, hibernation, migration, and molting, thus providing organisms with adaptive advantages. Although several hypotheses have been proposed, the genes that regulate circannual rhythms and the underlying mechanisms controlling long-term circannual clocks remain unknown in any organism. Here, we show a transcriptional program underlying the circannual clock in medaka fish (Oryzias latipes). We monitored the seasonal reproductive rhythms of medaka kept under natural outdoor conditions for 2 y. Linear regression analysis suggested that seasonal changes in reproductive activity were predominantly determined by an endogenous program. Medaka hypothalamic and pituitary transcriptomes were obtained monthly over 2 y and daily on all equinoxes and solstices. Analysis identified 3,341 seasonally oscillating genes and 1,381 daily oscillating genes. We then examined the existence of circannual rhythms in medaka via maintaining them under constant photoperiodic conditions. Medaka exhibited approximately 6-mo free-running circannual rhythms under constant conditions, and monthly transcriptomes under constant conditions identified 518 circannual genes. Gene ontology analysis of circannual genes highlighted the enrichment of genes related to cell proliferation and differentiation. Altogether, our findings support the "histogenesis hypothesis" that postulates the involvement of tissue remodeling in circannual time-keeping.

Regulation of PGC-1α of the Mitochondrial Energy Metabolism Pathway in the Gills of Indian Medaka (Oryzias dancena) under Hypothermal Stress.

Ectothermic fish exposure to hypothermal stress requires adjusting their metabolic molecular machinery, which was investigated using Indian medaka (Oryzias dancena; 10 weeks old, 2.5 ± 0.5 cm) cultured in fresh water (FW) and seawater (SW; 35‱) at room temperature (28 ± 1 °C). The fish were fed twice a day, once in the morning and once in the evening, and the photoperiod was 12 h:12 h light: dark. In this study, we applied two hypothermal treatments to reveal the mechanisms of energy metabolism via pgc-1α regulation in the gills of Indian medaka; cold-stress (18 °C) and cold-tolerance (extreme cold; 15 °C). The branchial ATP content was significantly higher in the cold-stress group, but not in the cold-tolerance group. In FW- and SW-acclimated medaka, the expression of genes related to mitochondrial energy metabolism, including pgc-1α, prc, Nrf2, tfam, and nd5, was analyzed to illustrate differential responses of mitochondrial energy metabolism to cold-stress and cold-tolerance environments. When exposed to cold-stress, the relative mRNA expression of pgc-1α, prc, and Nrf2 increased from 2 h, whereas that of tfam and nd5 increased significantly from 168 h. When exposed to a cold-tolerant environment, prc was significantly upregulated at 2 h post-cooling in the FW and SW groups, and pgc-1α was significantly upregulated at 2 and 12 h post-cooling in the FW group, while tfam and nd5 were downregulated in both FW and SW fish. Hierarchical clustering revealed gene interactions in the cold-stress group, which promoted diverse mitochondrial energy adaptations, causing an increase in ATP production. However, the cold-tolerant group demonstrated limitations in enhancing ATP levels through mitochondrial regulation via the PGC-1α energy metabolism pathway. These findings suggest that ectothermic fish may develop varying degrees of thermal tolerance over time in response to climate change. This study provides insights into the complex ways in which fish adjust their metabolism when exposed to cold stress, contributing to our knowledge of how they adapt.

Influence of phthalate and non-phthalate plasticizers on reproductive endocrine system-related gene expression profiles in Japanese medaka (Oryzias latipes).

Plasticizers containing phthalates have the potential to alter endocrine function in vertebrates. While non-phthalate plasticizers were previously considered to be environmentally friendly and safe, our research team discovered that bis-(2-ethylhexyl) adipate (DEHA) and acetyl tributyl citrate (ATBC) disrupt thyroid hormones in Japanese medaka (Oryzias latipes). We assessed reproductive- and estrogen-responsive gene expression patterns in Japanese medaka to determine whether the phthalate plasticizers bis-(2-ethylhexyl) phthalate (DEHP, positive control) and the non-phthalate plasticizers DEHA and ATBC disrupt endocrine signaling. The results showed that the levels of choriogenin H (chgH) and vitellogenin (vtg) genes increased after exposure to DEHP and ATBC, suggesting that these plasticizers may have estrogenic activity. Exposure to DEHP and DEHA resulted in the upregulation of kisspeptin (kiss), gonadotropin-releasing hormone (gnrh), and follicle-stimulating hormone beta (fshß) genes, suggesting that these plasticizers may interfere with reproductive function. To the best of our knowledge, this is the first study to demonstrate that the non-phthalate plasticizers DEHA and ATBC can disrupt reproduction-related hormonal activity in fish.

Expression of secretory calcium-binding phosphoprotein (scpp) genes in medaka during the formation and replacement of pharyngeal teeth.

BACKGROUND: Analyses of tooth families and tooth-forming units in medaka with regard to tooth replacement cycles and the localization of odontogenic stem cell niches in the pharyngeal dentition clearly indicate that continuous tooth replacement is maintained. The secretory calcium-binding phosphoprotein (scpp) gene cluster is involved in the formation of mineralized tissues, such as dental and bone tissues, and the genes encoding multiple SCPPs are conserved in fish, amphibians, reptiles, and mammals. In the present study, we examined the expression patterns of several scpp genes in the pharyngeal teeth of medaka to elucidate their roles during tooth formation and replacement. METHODS: Himedaka (Japanese medaka, Oryzias latipes) of both sexes (body length: 28 to 33 mm) were used in this study. Real-time quantitative reverse transcription-polymerase chain reaction (PCR) (qPCR) data were evaluated using one-way analysis of variance for multi-group comparisons, and the significance of differences was determined by Tukey's comparison test. The expression of scpp genes was examined using in situ hybridization (ISH) with a digoxigenin-labeled, single-stranded antisense probe. RESULTS: qPCR results showed that several scpp genes were strongly expressed in pharyngeal tissues. ISH analysis revealed specific expression of scpp1, scpp5, and sparc in tooth germ, and scpp5 was continually expressed in the odontoblasts of teeth attached to pedicles, but not in the osteoblasts of pedicles. In addition, many scpp genes were expressed in inner dental epithelium (ide), but not in odontoblasts, and scpp2 consistently showed epithelial-specific expression in the functional teeth. Taken together, these data indicate that specific expression of scpp2 and scpp5 may play a critical role in pharyngeal tooth formation in medaka. CONCLUSION: We characterized changes in the expression patterns of scpp genes in medaka during the formation and replacement of pharyngeal teeth.

A gene regulatory network combining Pax3/7, Sox10 and Mitf generates diverse pigment cell types in medaka and zebrafish.

Neural crest cells generate numerous derivatives, including pigment cells, and are a model for studying how fate specification from multipotent progenitors is controlled. In mammals, the core gene regulatory network for melanocytes (their only pigment cell type) contains three transcription factors, Sox10, Pax3 and Mitf, with the latter considered a master regulator of melanocyte development. In teleosts, which have three to four pigment cell types (melanophores, iridophores and xanthophores, plus leucophores e.g. in medaka), gene regulatory networks governing fate specification are poorly understood, although Mitf function is considered conserved. Here, we show that the regulatory relationships between Sox10, Pax3 and Mitf are conserved in zebrafish, but the role for Mitf is more complex than previously emphasized, affecting xanthophore development too. Similarly, medaka Mitf is necessary for melanophore, xanthophore and leucophore formation. Furthermore, expression patterns and mutant phenotypes of pax3 and pax7 suggest that Pax3 and Pax7 act sequentially, activating mitf expression. Pax7 modulates Mitf function, driving co-expressing cells to differentiate as xanthophores and leucophores rather than melanophores. We propose that pigment cell fate specification should be considered to result from the combinatorial activity of Mitf with other transcription factors.

Wnt4a Is Indispensable for Genital Duct Elongation but Not for Gonadal Sex Differentiation in the Medaka, Oryzias latipes.

In most vertebrates, the oviducts and sperm ducts are derived from the Müllerian ducts and Wolffian ducts, respectively. However, in teleosts, the genital ducts are formed by the posterior extension of gonads in both sexes. Whether the genital ducts of teleosts are newly evolved organs or variants of Müllerian ducts is an important question for understanding evolutionary mechanisms of morphogenesis. One of the genes essential for Müllerian duct formation in mice is Wnt4, which is expressed in the mesenchyme and induces invagination of the coelomic epithelium and its posterior elongation. Here, we addressed the above question by examining genital duct development in mutants of two Wnt4 genes in the medaka (wnt4a is orthologous to mouse Wnt4, and wnt4b is paralogous). The wnt4b mutants had a short body but were fertile with normal genital ducts. In contrast, both male and female wnt4a mutants had their posterior elongation of the gonads stopped within or just outside the coelom. The mutants retained the posterior parts of ovarian cavities or sperm duct primordia, which are potential target tissues of Wnt4a. The gonads of female scl mutants (unable to synthesize sex steroids) lacked these tissues and did not develop genital ducts. Medaka wnt4a was expressed in the mesenchyme ventral to the genital ducts in both sexes. Taken together, the data strongly suggest that the mouse Müllerian ducts and the medaka genital ducts share homologous developmental processes. Additionally, the wnt4a or wnt4b single mutants and the double mutants did not show sex-reversal, implying that both genes are dispensable for gonadal sex differentiation in the medaka.

High-speed system to generate congenic strains in medaka.

The congenic strain, an inbred strain containing a small genomic region from another strain, is a powerful tool to assess the phenotypic effect of polymorphisms and/or mutations in the substituted genomic region. Recent substantial progress in the genetic studies of complex traits increases the necessity of congenic strains and, therefore, a quick breeding system for congenic strains has become increasingly important in model organisms such as mouse and medaka. Traditionally, more than ten generations are necessary to produce a congenic strain. In contrast, a quick method has been reported previously for the mouse, in which the use of genetic markers reduces the required number of backcross generations to about a half that of the traditional method, so that it would take around six generations to obtain a congenic strain. Here, we present an even quicker congenic production system, which takes only about four generations. The system can produce medaka congenic strains having part of the HNI-II (an inbred medaka strain derived from the northern Japanese population, Oryzias sakaizumii) genome in the HdrR-II1 (another inbred strain from the southern Japanese population, O. latipes) background. In this system, the availability of frozen sperm and genotype data of the BC1 male population makes it possible to start marker-assisted congenic production after obtaining the BC2 population. Our evaluation revealed that the system could work well to increase the percentage of recipient genome as expected, so that a congenic strain may be obtained in about one year.

Altered Visual Function in Short-Wave-Sensitive 1 (sws1) Gene Knockout Japanese Medaka (Oryzias latipes) Larvae.

Visual perception plays a crucial role in foraging, avoiding predators, mate selection, and communication. The regulation of color vision is largely dependent on opsin, which is the first step in the formation of the visual transduction cascade in photoreceptor cells. Short-wave-sensitive 1 (sws1) is a visual pigment that mediates short-wavelength light transduction in vertebrates. The depletion of sws1 resulted in increased M-opsin in mice. However, there is still no report on the visual function of sws1 in teleost fish. Here, we constructed the sws1 knockout medaka using CRISPR/Cas9 technology. The 6 dph (days post-hatching) medaka sws1-/- larvae exhibited significantly decreased food intake and total length at the first feeding stage, and the mRNA levels of orexigenic genes (npy and agrp) were significantly upregulated after feeding. The swimming speed was significantly reduced during the period of dark-light transition stimulation in the sws1-mutant larvae. Histological analysis showed that the thickness of the lens was reduced, whereas the thickness of the ganglion cell layer (GCL) was significantly increased in sws1-/- medaka larvae. Additionally, the deletion of sws1 decreased the mRNA levels of genes involved in phototransduction (gnb3b, grk7a, grk7b, and pde6c). We also observed increased retinal cell apoptosis and oxidative stress in sws1 knockout medaka larvae. Collectively, these results suggest that sws1 deficiency in medaka larvae may impair visual function and cause retinal cell apoptosis, which is associated with the downregulation of photoconduction expression and oxidative stress.

Characterization of Fibrodysplasia Ossificans Progessiva relevant Acvr1/Acvr2 Activin receptors in medaka (Oryzias latipes).

Activin and Bone Morphogenetic Protein (BMP) signaling plays crucial roles in vertebrate organ formation, including osteo- and angiogenesis, and tissue homeostasis, such as neuronal maintenance. Activin and BMP signaling needs to be precisely controlled by restricted expression of shared receptors, stoichiometric composition of receptor-complexes and presence of regulatory proteins. A R206H mutation in the human (hs) BMP type I receptor hsACVR1, on the other hand, leads to excessive phosphorylation of Sons of mothers against decapentaplegic (SMAD) 1/5/8. This in turn causes increased inflammation and heterotopic ossification in soft tissues of patients suffering from Fibrodysplasia Ossificans Progressiva (FOP). Several animal models have been established to understand the spontaneous and progressive nature of FOP, but often have inherent limitations. The Japanese medaka (Oryzias latipes, ola) has recently emerged as popular model for bone research. To assess whether medaka is suitable as a potential FOP animal model, we determined the expression of Activin receptor type I (ACVR1) orthologs olaAcvr1 and olaAcvr1l with that of Activin type II receptors olaAcvr2ab, olaAcvr2ba and olaAcvr2bb in embryonic and adult medaka tissues by in situ hybridization. Further, we showed that Activin A binding properties are conserved in olaAcvr2, as are the mechanistic features in the GS-Box of both olaAcvr1 and olaAcvr1l. This consequently leads to FOP-typical elevated SMAD signaling when the medaka type I receptors carry the R206H equivalent FOP mutation. Together, this study therefore provides experimental groundwork needed to establish a unique medaka model to investigate mechanisms underlying FOP.

Environmental RNA as a Noninvasive Tool for Assessing Toxic Effects in Fish: A Proof-of-concept Study Using Japanese Medaka Exposed to Pyrene.

Although environmental RNA (eRNA) is emerging as a noninvasive tool to assess the health status of aquatic macroorganisms, the potential of eRNA in assessing chemical hazards remain largely untested. In this study, we investigated the ability of eRNA to detect changes in gene expression in Japanese medaka fish (Oryzias latipes) in response to sublethal pyrene exposure, as a model toxic chemical. We performed standardized acute toxicity tests and collected eRNA from tank water and RNA from fish tissue after 96 h of exposure. Our results showed that over 1000 genes were detected in eRNA and the sequenced read counts of these genes correlated with those in fish tissue (r = 0.50). Moreover, eRNA detected 86 differentially expressed genes in response to pyrene, some of which were shared by fish RNA, including the suppression of collagen fiber genes. These results suggest that eRNA has the potential to detect changes in gene expression in fish in response to environmental stressors without the need for sacrificing or causing pain to fish. However, we also found that the majority of sequenced reads of eRNA (>99%) were not mapped to the reference medaka genome and they originated from bacteria and fungi, resulting in low sequencing depth. In addition, eRNA, in particular nuclear genes, was highly degraded with a median transcript integrity number (TIN) of <20. These limitations highlight the need for future studies to improve the analytical methods of eRNA application.

High-fat diet in early life triggers both reversible and persistent epigenetic changes in the medaka fish (Oryzias latipes).

BACKGROUND: The nutritional status during early life can have enduring effects on an animal's metabolism, although the mechanisms underlying these long-term effects are still unclear. Epigenetic modifications are considered a prime candidate mechanism for encoding early-life nutritional memories during this critical developmental period. However, the extent to which these epigenetic changes occur and persist over time remains uncertain, in part due to challenges associated with directly stimulating the fetus with specific nutrients in viviparous mammalian systems. RESULTS: In this study, we used medaka as an oviparous vertebrate model to establish an early-life high-fat diet (HFD) model. Larvae were fed with HFD from the hatching stages (one week after fertilization) for six weeks, followed by normal chow (NC) for eight weeks until the adult stage. We examined the changes in the transcriptomic and epigenetic state of the liver over this period. We found that HFD induces simple liver steatosis, accompanied by drastic changes in the hepatic transcriptome, chromatin accessibility, and histone modifications, especially in metabolic genes. These changes were largely reversed after the long-term NC, demonstrating the high plasticity of the epigenetic state in hepatocytes. However, we found a certain number of genomic loci showing non-reversible epigenetic changes, especially around genes related to cell signaling, liver fibrosis, and hepatocellular carcinoma, implying persistent changes in the cellular state of the liver triggered by early-life HFD feeding. CONCLUSION: In summary, our data show that early-life HFD feeding triggers both reversible and persistent epigenetic changes in medaka hepatocytes. Our data provide novel insights into the epigenetic mechanism of nutritional programming and a comprehensive atlas of the long-term epigenetic state in an early-life HFD model of non-mammalian vertebrates.